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. Abstract 

o 

We apply Feynman's principle, "The same equations have the same solutions", to Kepler's problem 
and show that Newton's dynamics in a properly curved 3-D space is identical with that described by 
Einstein's theory in the 3-D optical geometry of Schwarzschild's spacetime. For this reason, rather 
unexpectedly, Newton's formulae for Kepler's problem, in the case of nearly circular motion in a 
static, spherically spherical gravitational potential accurately describe strong field general relativistic 
effects, in particular vanishing of the radial epicyclic frequency at r = r ms . 
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^ ■ 1. Introduction: observed QPOs as a test of strong gravity 



The standard accretion disk theory (Shakura and Sunyaev, 1973) assumes that matter in accretion 
disks around black holes and neutron stars moves on nearly circular, nearly Keplerian (i.e., nearly 
geodesic) orbits. Theory predicts several strong field general relativistic effects that should follow 
directly from this assumption, but none of them has so far been clearly detected. Most recently, 
we have found a new effect of this type (Kluzniak and Abramowicz, 2000, 2001, 2002). It concerns 
the QPOs, quasi periodic oscillations with kilohertz frequencies, observed as variations in the X-ray 
luminosity of accreting neutron stars and black holes. QPOs often occur in coupled pairs manifested 
as characteristic double peaks in the variability power spectra. According to us, observed frequencies 
at the double peaked QPOs are directly related to orbital epicyclic frequencies, vertical and radial, 
which are in a 3 : 2 parametric resonance. The resonance occurs in super strong gravity, just several 
gravitational radii outside the central black hole or neutron star, at a precisely determined radius. 
Knowing that, one may sharply constrain the global parameters of the source, as it was done, for 
example, for the Kerr angular momentum parameter in the black hole "candidate" GRO J 1655-40 
(Abramowicz and Kluzniak, 2001) that shows two QPOs at 450 Hz and 300 Hz (Strohmayer, 2001). 

In strong Einstein's gravity the radial epicyclic frequency £l r is smaller than the Keplerian orbital 
frequency Qk, and at the radius of the marginally stable orbit Q r = 0. In weak Newton's gravity 
one has Q r = fix ^ 0- These very different behaviours are attributed by many authors to a "non- 
linearity" of Einstein's gravity. We discuss here a simpler and more proper explanation: the difference 
is due only to the curvature of the three dimensional space. 
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Our explanation consists of three steps. First, we recall the relevant Einstein's equations. We write 
them in a particular form, consistent with the optical geometry of space. Second, we re-derive well- 
known, standard Newton's equations using a particular notation, and prove that they are identical 
in form with the corresponding Einstein's equations in optical geometry. Third, we apply Feynman's 
principle, "The same equations have the same solutions" (e.g. Feynman et al., 1989), and derive the 
formula for the epicyclic radial frequency that is valid in both Einstein's and Newton's gravity. 

2. Optical geometry in Schwarzschild spacetime 

The general static, spherically symmetric metric can be written in a particular form, 

ds 2 = e 2 * {-c 2 dt 2 + [drl + f 2 (dd 2 + sin 2 tidip 2 )] } . (2.1) 
In the specific case of Schwarzschild geometry one has, 

e 2* = l _r_G du = e -2* d ^ ~2 =e -2* r 2^ (22) 

r 

where tq = 2GM/c 2 is the gravitational radius of the central body with the mass M. The 3-D metric 
of optical geometry was introduced by Abramowicz, Carter and Lasota, (1988). It corresponds to the 
part of (2.1) in square brackets, 

^optical = dr 2 , + f 2 ((M 2 + sin 2 Mtp 2 ) , (2.3) 

In terms of optical geometry the Schwarzschild quantities (2.2) have a very clear geometrical meaning. 
Obviously, r* is the geodesic radius, and f is the circumferential radius of the nested, concentric spheres 
r = const that generate the spherical symmetry of the metrics (2.1) and (2.3). 

Let us introduce the covariant derivative operator V 1 in the optical geometry (2.3). From the Ru = 
component of Einstein's field equations one easily derives, 

V 2 $ = V l V,$ = 0, $ = C - (e 2 * - 1) (2.4) 

This optical geometry equation for $ is linear and identical with Laplace's equation that gravitational 
potential obeys in Newton's theory. This, together with the asymptotic behaviour of $ for large r, 
namely $ = —GM/r, suggests that <!> should be called the gravitational potential. We shall see later 
more reasons to do so. 

Equation of geodesic motion along a great circle on a r = const sphere, and with a constant speed v = 
c/?7, takes the form, valid both in the optical geometry, and in the full 4-D spacetime (Abramowicz, 
Carter and Lasota, 1988), 

v 2 

a i = C 2 V ^ X i = / 2 M 

n v ; 

Here a 1 is the four-acceleration, 7 = (1 — /? 2 ) -1 / 2 is the gamma Lorentz factor, A* = VV* is the 
first normal to the circle (orthogonal vector, of unit length in the metric of eq. [2.3]), and 7Z is the 
curvature radius of the sphere. (For simplicity, one may consider a particular great circle, located at 
the equatorial plane, $ = tt/2, but our arguments are valid in a general case). 
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The conserved angular momentum equals (Abramowicz, Carter and Lasota, 1988), 



C = we*. 



(2.6) 



We use (2.6) to write the final formula, 



C 2 ■ 

V l $ A* 

TZr 



= 0. 



(2.7) 



In these calculations one may use a convenient relation, 



dr 
dr* 



r 

n' 



(2.8) 



For light trajectories v = oo, and from (2.5) it follows that also 1Z = oo, which means that light may 
go round a circular trajectory if and only if this trajectory is a geodesic circle in optical geometry — 
a conclusion that follows also from Fermat's principle: light rays coincide with geodesic trajectories in 
optical geometry (Abramowicz, Carter and Lasota, 1988; Abramowicz, 1994). 

Equations (2.4) and (2.7) that we have recalled from Einstein's theory will be compared with corre- 
sponding Newton's equation that we derive next. 

3. Newton's dynamics in a curved 3-D space 

Newton's dynamics is usually considered in 3-D Euclidean space, but its generalization to a curved 
space is trivial. Indeed, the only issue that is important in the present context is a careful distinction 
between the three radii of a sphere: geodesic radius r*, circumferential radius f, and curvature radius 
1Z. In Newton's theory, with Euclidean geometry assumed, one has r* = f = 1Z, but assuming 
Euclidean geometry is not necessary in Newton's dynamics, and one could easily distinguish the three 
different radii in all calculations. 

Indeed, let us consider a curve in space defined by, 



x % = x z (s), 



(3.1) 



where x % are coordinates in a coordinate system, s is the length along the curve, and Latin indices 
run through 1, 2, 3. The velocity v l is defined by, 



dx 1 



(3.2) 



where t is the absolute time. From this definition it follows that 



dx 1 dx 1 ds 
dt ds dt 



(3.3) 



where t % = dx l /ds is a unit tangent vector to the curve (3.1), and v = ds/dt is the speed along the 
curve. 

Acceleration is defined as 



3 



and from this definition it follows 



a* = — = -r(r l w) = -r-r(T l «) = t> 2 — + ur l — = - — A + t 1 — — , (3.5) 
dt dV ' dtds y 1 ds ds K ds v 7 

where —A 1 is the first normal to the curve (3.1), 1Z is the curvature radius of the curve (3.1) and Sk 
is the kinetic energy per unit mass. Because £k is obviously constant for a circular motion with a 
constant speed, we may write that for such motion, 

v 2 ■ 

a 1 = -— A\ (3.6) 
Newtonian dynamics is based on the second law, 

F l = ma\ (3.7) 

where F l is the applied force and m the mass. Because we are interested here in a "Keplerian" motion, 
with 

F % = -mV l $ (3.8) 

being the gravitational force (equal to the gradient V* of gravitational potential $) the second law 
takes the form, 

v 2 • 

V l $ = —\\ (3.9) 
Using Newton's formula for angular momentum, 

C = vr, (3.10), 

we write finally, 

C 2 ■ 

V l $ = — A*. (3.11) 

This is identical with Einstein's equation (2.7). 

The gravitational potential obeys Laplace's equation, 

V l V,$ = 0, (3.12) 

which is identical with Einstein's equation (2.4). We thus have completed the second step, showing 
that for circular motion in spherical potential, Einstein's equations in optical geometry, and Newton's 
equations are the same. 
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4. The same equations have the same solutions 



Einstein's and Newton's equations for Kepler's circular motion in a spherical potential have the same 
form, 

v^ = ^y, (4.1) 

V l Vi$ = 0, (4.2) 

and the physical and geometrical meaning of all quantities appearing in them is the same. 

Let us integrate Laplace's equation (4.2) in the volume V between two equipotential surfaces, <Si, 
defined by <E>i = const, and S2, defined by $2 = const, and use the Gauss theorem, 

/ V\V^)dV= [ (V^)XldS= [ (VityX'dS- [ (V i $)\ i dS = 0, (4.3) 
Jv Js 1 +s 2 Js ± Js 2 

where A* is the vector orthogonal to the surface S — obviously, A* = A 1 at <Si and A* = —A 1 at 1S2. 
Let us fix a position of <Si , and denote 

f (VityX'dS = C. (4.4) 

If we now change the position of £2 then, because (4.3) holds independently of the position of 1S2, and 
(Vi<5)A l is constant on each equipotential surface, one must conclude that 



(Vi$)A l J dS = C = const, 



(4.5) 



for any equipotential surface S. Taking into account that r is the circumferential radius of the sphere, 
and therefore 



J dS = 4irf 2 , (4.6) 



one concludes that 



with C = AttGM following from the asymptotic behaviour at f — > 00. 
Thus, finally, the second law takes the form, 

GM C 2 



f 2 r 2 1Z J 

and from this we derive the formula for the Keplerian angular momentum distribution, 
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(4.8) 



C 2 = GMK. 



(4.9) 



The above formula allows a novel and interesting interpretation of the Keplerian angular momentum, 
as the geometrical mean of the gravitational radius of the gravitating center tq = 2GM/c 2 , and of 
the curvature radius 1Z of the particle trajectory, 



£ = c^2r G K. (4.10) 

One knows that the circular photon trajectory is located at r = (3/2)rc- This means that 1Z = oo 
there. But 1Z = oo also for r = oo. This means, that somewhere in the range [(3/2)rc oo] the 
curvature radius 7Z must have a minimum. Because C 2 = GM1Z, the angular momentum C also has 
a minimum at the same radius. It must be the radius r = r ms = 3rc of the mariginally stable orbit. 

To see this in a more quantitative way, let us now consider a small radial perturbation of Einstein's 
equation (2.7), which is equivalent to (4.1). We assume that a particle with a fixed angular momentum 
C is displaced from its original orbit at r* to a new one, at r* + <5r*. The perturbation introduces an 
unbalanced force and radial acceleration, 

d 2 d 2 1 dC 2 

e 2 *A-a* -e 4 *— (5r - e 4 *4 2 — Sr - — — 8r (411) 

6 Ka ~ e ds 2 * dt 2 * *' PTZdn ' [ ' 

where A denotes the total redshift factor, connected to the previously introduced Lorentz gamma 
factor by 

A 2 = 7 2 e" 2 * = (l - 1 . (4.12) 

From (4.11) we derive, finally, 

( dTZ f 2 \ 

4r - + *(*-^) & - = - (4 - 13) 

Here, each dot denotes a differentiation with respect to time of the observer at infinity, and 



/ GM\ 



1/2 



n K = Ce-^f-'A- 1 = (4.14) 
is the Keplerian orbital frequency, also observed at infinity. 

It is quite remarkable that equation (4.13) has identical form in Einstein's and Newton's theories — 
the well-known equation for a harmonic oscillator. Its eigenfrequency, 

n HSS)* (415) 

is obviously equal to the epicyclic frequency of small radial oscillations. When Q 2 > the radial 
epicyclic oscillations are stable, and when Q 2 < 0, they are unstable. Thus, mariginal stability occurs 
at the radius where dC/dr = = dlZ/dr which, in Schwarzschild geometry is at r = r ms . 
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In Newton's case the geometry is flat (Euclidean) and the geodesic, circumferential and curvature 
radii are equal, 



= f = n. 



(4.16) 



From these expression and (4.13) one derives, 



n 2 r = nl 



K 



in Newton's theory. 



(4.17) 



This means that the radial epicyclic frequency equals to the Keplerian orbital frequency, which is why 
Newton's orbits are closed ellipses. 

In Einstein's gravity, the geometry is curved and the geodesic, circumferential and curvature radii are 
all different, 



From these expressions and (4.13) one recovers the well-known formula, to our best knowledge first 
derived by Kato and Fukue (1980), 



This completes our point: Einstein's and Newton's formulae (4.1) and (4.2) that describe dynamics in 
static, spherically symmetric gravity are the same in both theories (and in both theories linear). For 
this reason, the formula for radial epicyclic motion (4.13), derived directly from these equations, is also 
the same in both theories. However, this formula depends on geodesic, circumferential and curvature 
radii of circular trajectories. In the flat, Euclidean, geometry of Newton's space these radii are equal, 
and in the curved geometry of Einstein's space they are not. This geometrical difference alone, and 
not an often mentioned (but not present in this case) "non-linearity" of Einstein's equations, is the 
reason for the distinctively different Newton's and Einstein's predictions for the physical behaviour 
of small radial oscillations around Keplerian circular orbits in static spherically symmetric gravity. 

The same conclusion was previously reached by Abramowicz, Lanza, Miller and Sonego (1997) for 
weak gravity, who considered the perihelion of Mercury advance according to Newton's gravity in 
a properly curved space. It is rather surprising that the conclusion holds also for arbitrarily strong 
static spherically symmetric gravity. 
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